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I. INTRODUCTION 



In a standard perturbative QCD approach to semi-inclusive processes and in particular to Semi-Inclusive 
Deep Inelastic Scattering (SIDIS), factorization theorems [l| allow to extract soft hadronic wave functions 
from high energy reactions data. Such non-perturbative process-independent distributions obey QCD renor- 
malization group equations Q. In presence of a hard scale, set by the virtuality of the exchanged boson in a 
Deep Inelastic event, standard parton and fragmentation distributions predict, together with the corresponding 
process-dependent coefficient functions [!, H|, the semi-inclusive cross-sections. These distributions, basic in- 
gredients in almost nowaday QCD-calculations, are well suited for studying full inclusive process, such as Deep 
Inelastic lepton-hadron Scattering or Drell-Yan process in hadronic collisions. In the recent past however it has 
became increasingly clear the less inclusive distributions, either space-like or time-like, are necessary to deal 
with a variety of semi-inclusive processes. In particular fc t -dependent distributions acquired particular relevance 
and a great activity has been registered recently in this research field In the SIDIS case, for istance, final 
state hadrons are expected to have a sizeable transverse momentum due to both intrinsic motion of partons into 
hadrons Q and to the radiative process off the struck parton line [3, Q . Unfortunately transverse momentum 
is usually integrated over, loosing part of the information which is contained in the experimental cross-sections. 
For these reasons it would be highly desiderable to have the evolution equations for these extended fct-dependent 
distributions. Such evolution equations were first proposed in the time- like case in Ref. [1, Q and then recently 
extended in the space-like domain in Ref. [10| . In order to have a complete description of the semi-inclusive 
cross-sections in terms of the transverse momentum, such a generalization was also performed in the target 
fragmentation region by introducing properly modified [ic| fracture functions [TlT |. The basic idea behind the 
fct-dependent evolution equations can be summarized as follows. Let us consider parton emissions off a active, 
space-like, parton line. In the collinear limit, at each branching, the generated transverse momentum is neg- 
ligible. In this limit however fc t -ordered diagrams can be shown to give leading logarithmic enhancements to 
the cross-sections. Since such contributions can be resummed by DGLAP evolution equations @, at the end 
of the radiative process, the interacting parton could possibly have an appreciable transverse momentum. As 
a result, fct-dependent evolution equations therefore depend, in addition to standard longitudinal momentum 
fraction, also on transverse degree of freedom. When solutions to the evolution equations are used to reproduce 
the SIDIS transverse momentum spectrum, the predictions smooth interpolate from small to large transverse 
momenta, this being a signature of well known DIS scaling violations in semi-inclusive process. 

The aim of this work is to offer a preliminar phenomenological study of fct-dependent evolution equations and 
to compare it with available hadron production data in DIS current fragmentation region. All the predictions 
are given by a handful of phenomenological assumptions. However such predictions are not the result of a fit 
to data, and thus strengthen our confidence in the general framework offered in Ref. [lCj . 
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II. TRANSVERSE MOMENTUM DEPENDENT EVOLUTION EQUATIONS 



Ordinary QCD evolution equations at leading logarithm accuracy (LLA) resum terms of the type 
a™ log n (Q 2 / fj.jp) originating from quasi-collinear partons emission configurations, where ijl 2 f represents the fac- 
torization scale. Leading contributions are obtained when the virtualities of the partons in the ladder are 
strongly ordered. At each branching, the emitting parton thus acquires a transverse momentum relative to its 
initial direction. The radiative transverse momentum can be taken into account through transverse-momentum- 
dependent evolution equations, which in the time-like case read [3J: 

n2 dV^(z h ,Q 2 ,p±.) a s (Q 2 ) f 1 du 2 f d 2 q ± 2 2\r,h( z h n i z h \ m 
dQ 2 = 2tt J — p ii\ u > a *(Q )) J —^—5[u(l-u)Q -q* x )V^y—,Q ,p±- — q±J. (1) 

Fragmentation functions X>f (z?,,, Q 2 ,p±) of eq. ([I]) give the probability to find, at a given scale Q 2 , a hadron h 
with longitudinal momentum fraction Zh and transverse momentum p±^ relative to the parent parton i. Pij(u) 
are the time-like splitting functions which, at least at LL accuracy, can be interpreted as the probabilities to find 
a parton of type i inside a parton of type j and are expressed as a power series of the strong running coupling, 
Pij{u) = J2 n =o a s {Q 2 )Pif > ( u )- The order n of the expansion of the splitting function matrix Pij(u) actually 
sets the accuracy of the evolution equations. The radiative transverse momentum square q\ at each branching 
satisfies the invariant mass constraint q\ = u (1 — u) Q 2 . The transverse arguments of T>\ (zh, Q 2 ,p±) on r.h.s. 
of eq. ([T]) are derived taking into account the Lorentz boost of transverse momenta from the emitted parton 
reference frame to the emitting parton one, see the left panel of Fig. (JTj) - 




FIG. 1: Boost of transverse momenta. Left panel: a time-like off-shell parton generated in a hard process, the grey 
blob, emits a daughter parton and acquires a transverse momentum q± relative to its intial direction. The small blob 
symbolizes the iteration of such emissions. Right panel: the analogue as before in the space-like case. 



The unintegrated distributions fulfil the normalization: 

J d 2 p ± V , ;{z h ,Q 2 ,p A _)=V<t{z h ,Q 2 ). (2) 

This property garantees that we can recover ordinary integrated distributions from unintegrated ones. The 
opposite statement however is not valid since eq. ([T]) contains new physical information. In analogy to the 
time-like case we consider now a initial state parton p in a incoming proton P which undergoes a hard collision, 
the reference frame being aligned along the incoming proton axis. We thus generalize eq. ([1]) to the space-like 
case fl(ij |: 

rt2 dfp( x B,Q 2 ,k ± ) a s (Q 2 ) f 1 du 2 f d 2 q± 2 2 , Tj /is 2 k± - q±\ 

Parton distribution functions J- p (xb,Q 2 ,k±) in eq. ([3]) give the probability to find, at a given scale Q 2 , a 
parton i with longitudinal momentum fraction xb and transverse momentum k± relative to the parent hadron, 
see the right panel of Fig. (Q}. The unintegrated distributions fulfil a condition analogous to the one in eq. J2|), 
i.e. : 

J d 2 k ± T P (x B ,Q 2 , k±) = r p {x B ,Q 2 ) ■ (4) 

We note that the inclusion of transverse momentum does not affect longitudinal degrees of freedom since partons 
always degrade their fractional momenta in the perturbative branching process. 

The approach can also be extended in the target fragmentation region of semi-inclusive D1S (Iol | by in- 
troducing a /c t -dependent version of fracture functions pd| . The corresponding evolution equations for 
M Ph (x, k±, z,p±,Q 2 ) can be obtained The factorization properties of these distributions however has not 
been proven yet, at variance with the current fragmentation case whose factorization in terms of fc t -dependent 
distributions has been proven in Ref. (I2I ]. 
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up quark gluon up(x,Q2,k^) 




FIG. 2: Space-like evolution. Left and middle panel: average transverse momentum < k\ > generated in the evolution 

of the up-quark and gluon for three different scales: Q% = 5 GeV 2 (-), Q 2 = 10 GeV 2 ( ) and Q 2 = 20 GeV 2 (-■) . 

Right panel: the transverse spectrum of the up quark at fixed xb for three different scales as before. The solid line is 
the gaussian initial condition. The evolved distributions show a l/(k 2 L y dependence. 



III. PHENOMENOLOGY IN THE CURRENT FRAGMENTATION REGION 

The fct-dependent evolution equations, eq. ((T|) and eq. ([21), are solved by means of a finite difference method 
in the (2n/ + l)-dimcnsional space of quarks, antiquarks and gluons. As appropriate for a leading logarithmic 
calculations, we set splitting functions to their lowest order expansion. In this preliminar analysis we simulate 
light flavours only while heavy flavours are accounted for in only as virtual contributions in the LL running 
coupling constant, a s (Q 2 ). Convolutions on transverse and longitudinal variables in eq. |T]) and eq. §3§ are 
numerically performed on a bidimensional (x,kj_) grid. To achieve a faster convergence and minimize the size 
of the grid, non-linear spacing both in x and in k\ have been adopted. At each Q 2 -itcration, the normalization 
conditions, eq. © and eq. are checked to reproduce ordinary longitudinal distributions within a given 
accuracy. As in the longitudinal case, fcj-dependent distributions at a scale Q 2 > Qq are calculable if one 
provides a non-perturbative input density at some arbitrary scale Qq. In the following we assume the simplest, 
physically motivated ansatz, i.e. a longitudinal parton distribution function Fp(xB,Qo) \)M or fragmentation 
functions D^(zh,Qo) (3 times a gaussian transverse factor, motivated by the Fermi motion of partons in 
hadrons @: 

, 2 2 
- fc J_ ~P± 

Fp(x B , Qo, fcj_) = Fj,(x B , Qq) 6 . 2 ° - , Pfr(z h , Ql, k x ) = D^(z h , Q 2 ) 6 ^ i = q,q,g- (5) 

Before comparing to data, we would like to draw some general properties of the evolution and discuss the 
stiffness of the initial conditions, eq. ([5]). We focus on the space- like case and set the width < k\ > to a 
testing value of 0.25 GeV 2 both for quarks and gluons. The evolution then is performed from the initial scale 
Qq — 5 GeV^ 2 to Q 2 — 20 GeV 2 , see Fig. [2] In order to reduce the number of parameters, we assume a flavour- 
independent value for the average transverse momentum < k 2 ^ > . Such hypothesis is indeed too crude in the 
quark valence region. Furthcmore, as it appears in Fig. the evolution generates a a^-dependent amount of 
averaged transverse momentum, behaving like 

<fcl>=</c 2 !0 > x%, 7 <0, (6) 

even starting from a SB-independent distribution, eq. |5]). This behaviour is expected since the arguments of 
/ct-dependent distributions in the right hand side of eq. ((T|) and eq. © mix, as a result of transverse boost, 
longitudinal and transverse degree of freedom. We have checked that the factorized form of eq. ([5]) is not 
preserved under evolution and deviation from a guassian form into broader fc^-distributions, especially for the 
gluon, are observed. In the rightmost panel of Fig. ^ is clearly visible how the evolution turns the the gaussian 
transverse factor at the initial scale into a inverse power-like distributions in k\ at the final scale. It is also 
visible in the same plot a de-population effect in the k 2 ^ ~ Q 2 region according to strong ordering recipe built-in 
the evolution equations. From above arguments and since the factorization scale Qq, at which we suppose eq. (O 
to be valid, is arbitrary we conclude that a more refined analysis could use initial condition with a a^-dependent 
transverse factor. We note also that the solutions do not show any growth of < fc 2 > in the large xb limit. In 
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FIG. 3: Left panel: average transverse momentum < P£ ± > versus zf, 100 < W < 340 GeV , Q > 5 GeV, against 
predictions (solid line). Right panel: average transverse momentum < P 2 ± > versus W 2 for 0.2 < z < 1.0, Q 2 > 5 
GeV 2 , against predictions (solid line). Data from Ref. f2C| . 



the soft limit the fc t -dependent evolution equations can be shown to diagonalize in impact parameter phase [l5j 
by a joint Fourier-Mellin transform 0. As a result soft gluon resummation technique can be applied to leading 
and next-to-leading logarithmic accuracy [l6l . [TtI Il8j . Attaining these limitations in mind we compare in the 
following the outcome of fct-dependent evolution equations with charged hadron production data in the DIS 
current fragmentation region. In this case we are supported by a factorization theorem and the semi-inclusive 
cross-sections can be shown to factorize in terms of fct-dependent distributions [12] • With leading logarithmic 
accuracy the cross-sections reads 



d 5 a 



dx B dQ 2 dz h dQ 2 d 2 P, 



h± 



e *J d 2 k ± d 2 p ± S^(z h k ± +p ± - P hl _)T P (x B ,Q 2 ,k ± ,)V^(z h ,Q 2 ,p ± ) , (7) 



while the soft factor, present in the original factorization formula in Ref. [jjj], is dropped for phcnomenological 
purposes. The standard SIDIS variables are defined as xb — Q 2 /{2P ■ q) and Zh = (P ■ Ph)/{P 4 q) where 
P, Ph, q are respectively the four momenta of the incoming proton, outgoing hadron and virtual boson. At 
lowest order, the process-dependent coefficient function is omitted and set to unity. Factorizations scales are 
set to fj/p = (j. 2 -) = Q 2 and large logarithmic ratios of the type log(/xf, D /Q 2 ) occuring in the perturbative 
calculations are moved in fc t -dependent distributions are then resummed by evolution equations. We compare 
our predictions with data of Refs. [HI, |2(j. These data sets are differential in the variable of interest and cover 
a broad kinematical region, where DGLAP dynamics is supposed to be valid. We note that the theoretical 
predictions which reproduce the data in the original Refs. [la H3| are based on QCD-calculations of Ref. [2l| . 
Light flavours average transverse momenta and factorization scale are then set to 

< >= 0.25 GeV 2 , < v\,D q ,Dn >= 0-20 GeV 2 , Q 2 = 5 GeV 2 , (8) 

for distribution and fragmentation functions respectively, according to Ref. [23]. The parameters in Ref. [22j] 
are obtained by a fitting procedure to the low-P^ ± differential cross-sections of Ref. [l9| using the same initial 
condition as given in eq. ([5]). We note that gluons in eq. |J7]) are absent since do not directly couple with 
the virtual boson but enter indirectly the cross-sections due to quark-gluon mixing in the evolution equations. 
Gluon widths are however essentially unknown and for this reason, in this preliminar analysis, we set them 
equal to light flavours parameters. In order to verify that this choice does not affect the presented results, we 
have checked that a 20% variation of gluon widths does not alter significantly the predictions in the kinematical 
region of Refs. 0, H3], the overall effect being a slightly slope variation of the large- Ph± tail in Fig. 
The role of gluon and its transverse spectrum is however of special interest especially in HERA and LHC 
kinematics, and thus certainly deserves a separated study. We require both the time-like and space-like evolved 
fct-dependent distributions to satisfy, both for quarks and gluons, the normalization condition, eq. ([2]) and ([4|), 
in the kinematical range of data da HI with an accuracy set to 10%: 



d 2 k^T t {x,k^,Q 1 



EMC 



K(x,Q 2 ) . (9) 



The accuracy however could by increased properly thickening the simulation grid. In Fig. ([3]) we show the average 
transverse momentum < P 2 ± > compared to the predictions of eq. ([7]) properly normalized to the relevant 



5 



0.1<z<0.2 



0.2<z<0.4 



0.4<z<1.0 



H 

CM £ 

H 
"D 


"D 

D 











> 

0) 


o 
o 

CM 
V 



V 

o 

in 



v 
o 
o 

CM 

CM 
> 
CO 


o 
w 
co 

v 






> 

CO 


o 
in 

CO 

A 



FIG. 4: The normalized semi-inclusive cross-sections in bins of z and W 2 for charged hadron production in the current 
fragmentation region. Data from Ref. [l9| : 



inclusive cross-sections. In the left panel a rise of < P 2 ± > with z 2 is observed. Essentially this dependence is 

_ > +z 2 < k\ > . The slope 

of the data is roughly reproduced. On the right panel of Fig. ^ the < P 2 j_ > results obtained from eq. ([7]) as 

> spectrum shows a clear logarithmic 



guided by the <5^ 2 ^-function in eq. ([7|) which leads to the expectation < P^ ± >=< p 

h± 



a function of W 2 = Q 2 {1 — xb)/xb are compared to data. The < P% d 

dependence on W 2 and the predictions far underestimate the measured average transverse momentum. As can 
be seen in Ref. [13], the measured dependence of < P, 2 ± > on Q 2 is very mild, while the one on xb is steeper 
and of a kind shown in Fig. The fej-dependent evolution equations take care of the former while probably 
only a ss-dependent correction in the transverse factor in eq. ([5]) could solve the latter. The charged hadron 
production P^-diffcrcntial cross-sections, properly integrated in the relevant z and W 2 bins and normalized to 
the inclusive total cross-sections, is shown in Fig. Q. The main effect of evolution equations is of modifying the 
sharp tail of the guassian distributions at Qq, see Ref. |22j |. into a broader transverse power-like distributions. 
At fixed W 2 , a progressive broadening of the spectrum, according to the left panel of Fig. ([3]), is observed. At 
fixed z instead, the predictions fall more distant from data as long as W 2 increases, according with Fig. (J3j> - 
The overall agreement looks however encouraging since we have not performed any fit to the data, apart from 
fixing the transverse widths as already discussed. At high W 2 , in the low-P^ part of the spectrum, deviations 
in slope beetwen data and predictions are visible, signalating again the inadeguacy of a x_B-mdependent width. 
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The underestimation of the transverse spectrum at high P?j_ indicates instead that large angle parton emissions 
from fixed order matrix element are needed. We conclude that both a more accurate choice of the initial 
condition and the inclusion of next-to-leading corrections will lead thus to a better agreement of the predicted 
cross-sections with data. 



IV. CONCLUSIONS 

In this work, by using the factorization theorem of Ref. [T^, the charged hadron production cross-sections 
in the current fragmentation region has been computed within leading logarithmic approximation by using 
fct-dependent evolution equations. The obtained /cj-dependent distributions, due to resummation of soft and 
collinear parton emissions, reproduce the high P^_i_ tail of tranverse spectra and smoothly interpolate from the 
low to the high P£ ± regime without using any matching procedure between the two regions. A reasonable 
description of the data is obtained by only using default width values as proposed in Ref. [12]. This validates 
our approach as proposed in Ref. [101 ]. The impact of the intial conditions, eq. (0, is investigated and arguments 
for a ^-dependent transverse factor are given, along with hints suggesting the need of higher order corrections. 
In this work we do not emphasize gluon dynamics. This subject however is a central one, especially for HERA 
and LHC kinematics [23|, \2M > an d we deserve it for a separate study. 

We wish to conclude by listing two possible promising applications of the presented formalism. The fct- 
dependent evolution equations could be tested in Drell-Yan pair production cross-sections differential in the 
transverse momentum of the lepton pair. The present formalism could find interesting applications to polarized 
reactions and could be particularly fruitful, for istance, in the case of transversity ditributions [25[. 
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